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Abstract. We study the change in free and internal energy due to the presence of a heavy quark anti-quark 
pair in a thermal heat bath in QCD with 2-flavors of staggered quarks at finite temperature. We discuss 
string breaking below as well as screening above the transition. Similarities and differences to the quenched 
case are discussed. 
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1 Introduction 

The study of the fundamental forces between quarks is 
an essential key to the understanding of QCD and the 
occurrence of different phases at high temperatures (T) 
and/or non-zero quark chemical potential (ix). The free 
energy of a static quark anti-quark pair 7j , separated by 
distance r, is an important tool to analyze the r and T 
dependence of the QCD forces. Similar to the free energies 
also the internal energies have been introduced [2]|3| and 
are expected to play an important role in the discussion 
of quarkonia binding properties QEJS]. The structure of 
these observables in the short and intermediate distance 
regime, rT < 1, is relevant for the discussion of in- medium 
modifications of heavy quark bound states |7ll51IMTnilTT] 
which are sensitive to thermal modifications of the heavy 
quark potential [12] • Up to quite recently (13114*] . however, 
most of these discussions concerned quenched QCD. 

Several qualitative differences, however, are to be ex- 
pected when taking the influence of dynamical fermions 
into account; the phase transition in QCD will appear as 
an crossover rather than a 'true' phase transition with re- 
lated singularities in thermodynamic observables. More- 
over, in contrast to a steadily increasing confinement in- 
teraction in quenched QCD, in full QCD the strong inter- 
action below T c will show a qualitative different behavior 
due to the possibility of string breaking. Thus it is quite 
important to extend our developed concepts for the analy- 
sis of free energies and internal energies in quenched QCD 
Plll5lllt>j to the complex case of QCD. This has recently 
been done for 2- and 3-flavor QCD [HlfTTlfT^. 



2 Quark anti-quark free energy 

We will discuss here lattice results for the quark anti-quark 
free and internal energies in 2-flavor QCD (Nf = 2) using 



an improved staggered fermion action with quark mass 
m/T = 0.4. Any further detail on this study can be found 
in [13117118) . For information on the improved actions 
used in these simulations see [TT?ll2()ll21j . While in earlier 
studies of the heavy quark free energy in quenched |22) 
and full QCD the color averaged operators were ana- 
lyzed we discuss here the free and also the internal energies 
in the color singlet channel. The static quark sources are 
described by the Polyakov loop, W(x) = YIt=i Uo{x,t) 
with Uq(x,t) S 5*7(3) being defined on the lattice link 
in time direction. The free energies in the color singlet 
channel is then given by [TTl24| 



-Fx(r)/T + C = L Tl{W ( x)W ^0)) , 



(1) 



where r = \x\ and C is a suitably chosen renormalization 
constant. The operator used here to calculate Fi(r,T) is 
not gauge invariant. Our calculations thus have been per- 
formed in Coulomb gauge "2"""] • 



2.1 Renormalized free energy 

As noted above, the quark anti-quark free energy |Tj"} and 
the expectation value of the Polyakov loop suffer from di- 
vergences and require renormalization. We follow here the 
conceptual approach recently suggested for the quark anti- 
quark free energies in quenched QCD [2] and apply it to 
our lattice results in full QCD. First experiences with this 
renormalization prescription in QCD are reported in [131 
114) . Renormalization prescriptions for the Polyakov loop 
expectation value are discussed in 28,29 . 

At distances much smaller than the inverse tempera- 
ture, r <C 1/T, the dominant scale is set by r and the run- 
ning coupling will be controlled by this scale, r -C 1/Aqcd 
26 , 1*"| ■ ln this limit the free energies are dominated by 
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Fig. 1. Results for the renormalized color singlet free energies 
calculated in 2-flavor QCD for 0.75 < T < 2.00. The solid 
line represents the heavy quark potential at zero temperature, 
V(r), given in |26| . At T = string breaking is expected at 
distances about 1.3 — 1.5 fm |221- The corresponding energies, 
V(1.3 - 1.5 fm) ~ 1100 - 1285 MeV, are shown by the band. 



Fig. 2. The Debye screening radius, m{T), and the scale 
r m ed(T) as function of temperature in physical units obtained 
from 2-flavor QCD lattice studies (filled symbols). The open 
symbols correspond to calculations in quenched QCD 121 1 
Further analysis of these data with respect to quarkonium dis- 
sociation temperatures can be found in |30| . 



one-gluon exchange, i.e. are given by the heavy quark po- 
tential, 

4 air) 

F t {r,T) = V{r) ~ -3-^. (2) 

We have neglected here any constant contributions to the 
free energy which, in particular, will dominate the large 
distance behavior of the free energy. Moreover, we already 
anticipated here the running of the coupling with the dom- 
inant scale r. Following [2] we can use this property to 
fix the constant C in by matching F\ (r, T) to the 
zero temperature heavy quark potential at small distance, 
F\(r <C 1/T,T) ~ V(r). Once the free energy is fixed at 
small distances also the large distance behavior is fixed as 
no additional divergences get introduced at finite temper- 
ature. 

Our lattice results for F% (r, T) calculated in 2-flavor 
QCD are summarized in Fig.^at several temperatures in 
the vicinity of the phase transition. Due to renormaliza- 
tion the free energies coincide with the heavy quark po- 
tential, V(r), at small distances (V(r) is specified in |17l 
I18U26| ). When going to larger distances, however, ther- 
mal modifications become important and the different ef- 
fects from color screening (T > T c ) and string breaking 
{T ^ T c) can be studied. 

2.2 Color screening and string breaking 

As the free energies shown in Fig. ^ rapidly change from 
the zero temperature like behavior at small distances to an 
almost constant behavior at large distances we introduce 
here scales that may characterize this qualitative change. 
For this purpose we introduced in [2] a scale r rnec i de- 
fined through Foo(T) = V(r med ), where F^T) is the 
plateau value which the free energy approaches at large 



distances, F^T) = ]im r -> 00 Fi(r,T). This scale charac- 
terizes the typical distance at which string breaking and 
color screening become relevant. Alternatively one can at 
high temperature characterize screening properties of the 
QCD plasma (T > T c ) in terms of the Debye screening 
length, ru(T), which is commonly defined as the inverse 
of the screening mass, ro(T) = l/m£>(T). We also esti- 
mated this scale by extracting the screening masses non- 
perturbatively from the exponential fall off that the free 
energies show above T c at large distances. Our results for 
r D{T > T c ) and r me d are shown in Fig. [2] as function of 
T/T c and are compared to findings in quenched QCD 2 , 
lib] . In general we find that both scales decrease with in- 
creasing temperature and r me d(T) > rjy(T). We note that 
both scales are expected to behave like 1/gT in the high 
temperature limit. On the other hand, when comparing 
both scales to the findings in quenched QCD no or only 
little differences in the temperature range 1.2 < T/T c < 2 
could be identified. However, at temperatures close to and 
below the transition, T < 1.2T C , differences become quite 
apparent. These differences signal the qualitative change 
in the phase transition when changing from quenched (I s *- 
order transition) to full QCD (crossover). In QCD r me d is 
finite also below T c due to string breaking and shows a 
rapid increase with decreasing temperatures 1 . In fact, the 
distance where the string is commonly expected to break 
at T = [27], r med {T = 0) ~ 1.3 - 1.5 fm, is almost 
approached already at T ~ 0.8T C . The quark anti-quark 
free energies in 2-flavor QCD thus are expected to show 
only little deviations from V{r) at T = at T > 0.8T C . 

This is also evident when analyzing the plateau val- 
ues, Foo(T), summarized in Fig. [3] as function of T/T c 
and again compared to recent findings in quenched [3] 
and 3-flavor QCD ^3] for temperatures in the vicinity 

In quenched QCD F\ (r, T < T c ) signals strict confinement 
and r m ed is infinite. 
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Fig. 3. The plateau value of the quark anti-quark free energy, 
Foo(T), calculated in 2-flavor QCD as function of T/T c at tem- 
peratures in the vicinity and below the phase transition. The 
dashed band show again the string breaking energies at T = 0, 
V(rbreak), with r bre ak — 1-3 - 1.5 fm - The open sym- 
bols (T > T c ) correspond to Foo(T) in quenched QCD (from 
Ref. 0) and the crosses to 3-flavor QCD studies (from (14) 't . 



and below T c . 2 At T ~ 0.8T C F^T) still is compatibel 
with the T = values, V(r = oo) ~ 1.1 GeV given in 
8,27 . It, however, rapidly drops in the vicinity of the 
transition to about half of this value, F^iTc) ~ 570 MeV. 
Although some flavor dependence can clearly be identified 
when comparing results in 2- and 3-flavor QCD the value 
at T c is almost the same. A similar value is also found in 
quenched QCD just above T c . At very high temperatures 
perturbation theory j^J suggests also negative values for 
Foo(T), i.e. 

Foo(T) * -~m D (T)a(T) ~ ~0(g 3 T) . (3) 

Such a behavior has indeed been observed in lattice studies 
of quenched QCD at high temperatures 0E1E21- 



3 Quark anti-quark internal energy and 
entropy 

It has been argued that the quark anti-quark free en- 
ergy, Fi(r, T) = U x {r : T) - TS^r, T), contains non-trivial 
entropy contributions, i.e. S± = S±(r,T), which, in par- 
ticular, could make the analysis of thermal modifications 
of the finite temperature potential complicated • We 
briefly discuss here the quark anti-quark internal energy 3 , 



2 The comparison of F^T) performed here clearly depends 
on the relative normalization of the corresponding heavy quark 
potentials used for renormalization. We fixed these over all 
contributions such that the Cornell parameterization of V(r) 
in all cases (Nf = 0, 2, 3) contains no constant contribution. 

3 As we calculate the internal energy and entropy using the 
renormalized free energies, Ui(r,T) and Si(r,T) are properly 
fixed by construction. However, the quark anti-quark internal 




Fig. 4. The asymptotic free energy, Foo (line), and the internal 
energy, U x (filled symbols), in 2-flavor QCD as function of 
T/T c . The dashed band shows again V(r = oo) ~ 1100 - 1285 

MeV 151271. 



Ui(r,T) = -T 2 d(Fi/T)/dT, and entropy, Si(r,T) = 
-dF 1 /dT. 

It is important to realize from © that although the 
entropy contribution at (infinite) large distances, Soo(T), 
will vanish in the perturbative high temperature limit, 
the quantity T5 , 00 (T) will increase and dominate the dif- 
ference between free and internal energy, i.e. 

UooiT) - Foo(T) = TS^T) ~ +0(g 3 T) . (4) 

In the limit of zero temperature, however, the observable 
TSoo(T) is supposed to vanish as Soo(T) is a dimension- 
less quantity. A qualitative change of the temperature 
dependence of the observable T5' 00 (T) when going from 
T = to high temperatures, T — > oo, is not obvious and if 
present will demonstrate the phase change from confine- 
ment to deconfinement. 

Our results for -Foo^) and Uoo(T) are summarized in 
Fig. Has function of T/T c . The observable, TS' 00 (T) = 
Uoo{T) — FoaiT), can easily be deduced when compar- 
ing the values for Uoo(T) with the values for Foo(T). It 
can clearly be seen that Uoo(T) and T5' 00 (T) exhibits a 
sharp peak at T c . The internal energy falls from a value 
of UoJT) w 4000 MeV at T c to about half of this value 
at a temperature of T/T c ~ 0.8. While the temperature 
dependence of Uoo(T) is quite strong in the vicinity of 
the transition and markedly different from the free en- 
ergy, both quantities show a similar T-dependence for 
1.5 > T/T c > 2. This is also expected at lower tem- 
peratures (T < 0.7T C ) were U^T) and F^T) are ex- 



energy and entropy could also be seperately renormalized non- 
pcrturbatively at small distances, i.e. U\(r 1/T, T) ~ V(r) 
and Si(r < 1/T, T) ~ 0. 

It is interesting to note here that in contrast to Fi (r, T) and 
Ui(r, T), the the entropy is free from any finite renormalization 
at T = which, in particular, could introduce flavor dependent 
over-all constant contributions to the finite temperature ener- 
gies. 



4 



Olaf Kaczmarek and Felix Zantow: Static quark anti-quark free and internal energy in 2-flavor QCD 




0.5 1 1.5 2 



Fig. 5. The color singlet quark anti-quark internal energy 
above T c . The solid line represents the T=0 heavy quark po- 
tential, V(r). 

pected to smoothly approach the same value at T = 0, 
V(r = oo ) — 1100 - 1285 MeV (from (HIEZI ) - 

Finally we show in Fig.JSJthe r-dependence of the inter- 
nal energies above T c . Indeed Ui(r, T) agrees with the zero 
temperature potential at small distances while at interme- 
diate distances U\ (r, T) has a steeper slope and clearly lies 
well above Fi(r,T) (compare with Fig.0. Further details 
on our analysis of the finite temperature energies, F\(r, T) 
and U\ (r, T) , and the entropy will be given |17U18| . 

4 Conclusions 

We have discussed thermal modifications of quark anti- 
quark free and internal energies in the color singlet channel 
in QCD with dynamical quarks (Nf = 2) and compared 
our findings to similar studies in quenched and 3-flavor 
QCD. Although our comparison with the quenched case 
is in parts only on a qualitative level, we already leads 
to useful information for the study of heavy quark bound 
states in the plasma phase. At temperatures well above 
T c , i.e. 1.2 < T/T c < 2, no or only little differences ap- 
pear between results obtained in quenched and 2-flavor 
QCD. This might suggest that using thermal parameters 
extracted from energies in quenched QCD as input for 
model calculations is a reasonable approximation. This 
also supports direct investigations of heavy quarkonia in 
quenched lattice QCD 33 34 35]. We note, however, that 
most of our 2-flavor QCD results differ from quenched 
calculations at temperatures in the vicinity and below the 
transition. These differences could make a discussion of 
possible signals for the quark gluon plasma production 
based on quenched QCD results complicated when tem- 
peratures close to the transition become important. 
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